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Abstract. We present a framework, compliant with the general canonical principle of 
statistical mechanics, to define measures on the set of pure Gaussian states of continuous 
variable systems. Within such a framework, we define two specific measures, referred to 
as 'micro-canonical' and 'canonical', and apply them to study systematically the statistical 
properties of the bipartite entanglement of n-mode pure Gaussian states (as quantified by the 
entropy of a subsystem). We rigorously prove the "concentration of measure" around a finite 
average, occurring for the entanglement in the thermodynamical limit in both the canonical and 
the micro-canonical approach. For finite n, we determine analytically the average and standard 
deviation of the entanglement (as quantified by the reduced purity) between one mode and all 
the other modes. Furthermore, we numerically investigate more general situations, clearly 
showing that the onset of the concentration of measure already occurs at relatively small n. 
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1. Typical entanglement in quantum information theory 

Due to the exponentially increasing complexity of the Hilbert spaces of multiple constituents, 
a complete theoretical characterisation of the entanglement of general quantum systems of 
many particles turns out to be a daunting task (TJ. A viable approach towards such a 
characterisation consists in focusing on the "typical", statistical properties of the quantum 
correlations of multipartite systems, when the states of the system are assumed to be 
distributed according to a particular 'measure'. This strategy is firstly aimed at simplifying 
the problem at hand by restricting attention on the typical (and thus, in a sense to be precisely 
specified in the following, "overwhelmingly likely") features of the entanglement of a system 
whose state is apt to be described by the chosen measures. Furthermore, this kind of analysis 
is able to shed light on the general properties of the entanglement of physical systems. 

For finite dimensional quantum systems, a natural, 'uniform' measure on pure states 
emerges from the "Haar" measure of the unitary group (i.e., from the left- and right-invariant 
measure under application of any unitary transformation), whose elements allow to retrieve 
any state when applied to another given starting pure state. On such grounds, a well defined 
typical entanglement of finite dimensional systems can be addressed and analysed. Original 
studies in this direction were undertaken well before the development of the formal theory of 
entanglement developed in quantum information science (TJ: In 1978, Lubkin considered 
the expected entropy of a subsystem when picking pure quantum states at random from 
the uniform measure||2). Let us recall that the von Neumann entropy S — —trphxp of a 
subsystem in state p properly quantifies, for globally pure states, the entanglement between 
the subsystem and the remainder of the system. Lubkin showed that one expects this quantity 
to be nearly maximal. Pagels and Lloyd arrived later at the same qualitative conclusions, 
following an independent line of thought [3|. Their work was expanded by Page, who 
conjectured an exact formula for the average entropy of a subsystem S m) „ EJ, reading 



for a quantum system of Hilbert space dimension mn in a random pure state, and a subsystem 
of dimension m < n. This relation was later proven by Foong and Kanno Q. 

This general line of enquiry was revisited and considerably extended in the setting 
of quantum information theory by Hay den, Leung and Winter in Ref. J6), where they 
extensively studied the 'concentration of measure' around the average of the entanglement 
probability distribution with increasing n. They also pointed out that this study may provide 
a way of simplifying the theory of entanglement which contains a plethora of locally 
inequivalent classes. In other words, as already mentioned, restricting statements to the 
"typical entanglement" allows one to ignore several unessential complications. Recent results 
on the physical interpretation of Page's conjecture can be found in Refs. [8), where it 
is proven that a circuit of elementary quantum gates on a quantum circuit is expected to 
maximally entangle the state to a fixed arbitrary accuracy, within a number of gates that grows 
only polynomially in the number of qubits of the register. 

A further simplification in the analysis of the entanglement can be achieved by 
considering the typical entanglement of 'particularly relevant' (according to the specific 
problem at hand) subsets of states. For example, it has been found that 'stabilizer states' 
(a countable set of states playing a central role in quantum error correcting codes) are also 
typically maximally entangled |J9] [TOl, similarly to the set of all states. Such investigations 
are interesting per se, as they unveil the potential and limitations hidden in the adoption of 




(1) 



k=n+l 
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restricted classes of states and, furthermore, provide us with a more detailed understanding of 
the entanglement properties of the total state space. 

Here we expand these considerations into the realm of continuous variables, 

i.e. of quantum systems described by pairs of canonically conjugated observables with 
continuous spectra. Such systems, ranging from motional degrees of freedom 

of particles in first quantisation to bosonic fields in second quantisation, are ubiquitous 
to all areas of quantum physics, being prominent in quantum optics (as they embody the light 
field in second quantisation), atomic physics (notably, in the description of atomic ensembles), 
quantum field theory (as they encompass any bosonic field), in addition to their crucial role in 
molecular and atomic physics. 

Quantum systems described by operators with continuous spectra live in infinite 
dimensional Hilbert spaces. Therefore, a first naive try in this direction could be to take the 
infinite dimensional limit of Eq. ([TJ. Page showed that Eq. (T]i implies that S m ^ n ~ In m — 
for 1 <C m < n. Under that restriction, and noting that the maximal entanglement is In to, 
we can then make the observation that the ratio of the entanglement average to the maximum 
tends to unity as n — ► oo, but that the two quantities both diverge logarithmically. From 
this perspective, the entanglement could be said to be typically infinite in the continuous 
variable setting. Even though mathematically reasonable (in the limit's sense), this statement 
is definitely questionable, physically and practically. In fact, in any practical situation, one 
will deal with a finite total energy or with finite "temperatures" (both quantities will be defined 
precisely in our treatment), whereas infinitely entangled states require an infinite energy to be 
created. In this paper, we shall restrict our attention on pure Gaussian states whose typical 
entanglement we shall study under two different measures, both inspired by arguments of 
thermodynamic al natures, but apt to describe different situations. The previously mentioned 
divergences, still potentially emerging in the Gaussian setting, will be tamed by introducing 
proper prescriptions which will generally affect the energy of the system, by imposing either 
a sharp upper bound or an exponentially decaying distribution of energies (in this respect, see 
also Ref. ifTTIl for a detailed discussion of the impact of similar constraints on entanglement 
measures). We will show that such prescriptions induce the occurrence of a finite typical 
entanglement in the limit of an infinite number of total constituents. Even for a finite number 
of total degrees of freedom and finite upper bound to the energy - entailing the existence of a 
finite maximal entaglement - the typical entanglement concentrates around a value which is 
well distant from the allowed maximum. 

Notice that Gaussian states are certainly the most prominent class of states not only in 
quantum information with continuous variables, but also, more broadly, in quantum optics, 
as they can be generated and manipulated with relative ease (even as highly entangled states 
lfT^[T3l[T4l ). can be used for the implementation of quantum communication and information 
protocols lfl5l and serve as a powerful testing ground for the theoretical characterisation of 
entanglement properties lfl6l . 

The micro-canonical measure, which we will apply here to the study of the typical 
entanglement, has been already employed in the analysis of the quantum teleportation of 
Gaussian states with generic second moments ATI . In particular, the micro-canonical average 
quantum fidelity and a corresponding "classical threshold", have been evaluated for the 
teleportation of states with null first moments and arbitrary second moments under the 
standard continuous variable teleportation protocol (see, e.g., lfl5ll for a description of the 
scheme). 

Let us finally mention that a definition of micro-canonical average entanglement has been 
very recently addressed for finite dimensional systems as well ifTHl . with a major emphasis on 
the possibility of reducing time-averages to ensemble-averages. 
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This paper is organised as follows. In Section [2] we review some preliminary facts 
about Gaussian states and set the notation. In Section [3] we review the definition of the 
micro-canonical measure on pure Gaussian states, already introduced in ifTTl , completing 
our previous analysis with the inclusion of comments and mathematical details previously 
omitted, and we extend the existing framework to encompass a 'canonical' measure as 
well. Section |4] contains a rigorous proof of the 'concentration of measure', common to the 
two measures introduced here, i.e. of the fact that the entanglement probability distribution 
concentrates in the thermodynamical limit, around a finite 'thermal' average, away from the 
allowed maximum. Even though this result had been anticipated in Ref. IfTTl . this proof is 
original and adds further insight into the matter. Section[5]presents a detailed study about the 
typical entanglement of pure Gaussian states with finite number of degrees of freedom, where 
both analytical findings and numerical evidences are reported. Conclusions and outlook are 
found in Section [6] Three appendices complement the work, one of which ( Appendix B I 
contains the most technical steps needed to prove the concentration of measure, while in 
Appendix A a specific Baker-Campbell-Hausdorff relation is derived, and in Appendix C a 



derivation of the expression of the maximal entanglement for given energy is presented. 



2. Preliminaries 



We consider bosonic continuous variable (CV) quantum mechanical systems described 
by n pairs of canonically conjugated operators {xj, pj} with continuous spectra, like 
motional degrees of freedom of particles in first quantisation or bosonic field operators 
in second quantisation. Grouping the canonical operators together in the vector R = 
(xi, . . . , x n ,pi, . . . ,p n ) T allows to express the canonical commutation relations (CCR) as 
[Rj,Rk] = 2i Cljk, where the 'symplectic form' f2 is defined as 




n and 1„ standing for the null and identity matrix in dimension n. 

Any state of an rt-mode CV system is described by a positive, trace-class operator q. For 
any state q, let us define the 2rc-dimensional vector of the expectation values ("first moments") 
of the canonical operators R (with entries Rj) as 

Rj = Tr [RjQ] 

and the 2n x 2n matrix of second moments, or "covariance matrix", cr (with entries aij) as 

a id = Tr [{R 1 ,R j }q}/2 - Tr [R lQ ]Tr [RjQ] . 

Also, throughout the paper, we will refer to the 'energy' of a state q as to the expectation 
value of the operator Ho = Sj=i(^j + Pj)- This definition corresponds to the energy of a 
free electromagnetic field in the optical scenario (and to decoupled oscillators in the general 
case). In our convention, as determined by the factor 2 appearing in the CCR, the vacuum of a 
single mode has covariance matrix equal to the identity (thus simplifying significantly several 
expressions), with energy 2 (the adopted energy unit is Hu>/4 for a mode of frequency uf). The 
energy is determined by first and second moments according to 

Tr(gHo) = Tr(<r) + ||i?|| 2 , 

where ||i?|| is the usual euclidean norm of the vector R. 
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Gaussian states are defined as the states with Gaussian characteristic functions and quasi- 
probability distributions, defined over a phase space analogous to that of classical Hamiltonian 
dynamics. As well known, a pure state \i/j)g is Gaussian if and only if it can be obtained by 
transforming the vacuum |0) under an operation generated by a polynomial of the second 
order in the canonical operators. In formulae (up to a negligible global phase factor): 

\4>)g = G A ,b\0) = e« ATAA+AT V |0) , (2) 
where A and b are, respectively, a real 2n x 2n matrix and a real 2n-dimensional vector. 
Because of the CCR and of the unitarity of Gam = c l ( R AR+R b \ A can be chosen 
symmetric, without loss of generality, while b is a generic real vector. A and b determine 
the CM and the second moments of the Gaussian state, thus completely determining it. 
The unitary operator can always be rewritten as (see Appendix A i 

G Ab = S rTa ^ S kTMb ~> , (3) 



for some matrix M (as shown in Appendix A M = CIA (l2n — e )/4 for invertible 



A's). First and second order operations can thus be generally 'decoupled'. 

First order operations correspond to local displacements in phase space. While such 
operations do not affect local entropies (and thus the entanglement) of multipartite states, they 
do affect the energy of the states, which will play a central role in what follows. Moreover, 
let us notice that the group of these transformations is non-compact, being isomorphic to 
the abelian R 2 ™ under the addition composition rule. In the following, we will show how 
the first moments can be consistently incorporated in the presented framework. However, 
because their inclusion in the study of the statistical properties of the entanglement is just a 
technicality (adding no significant insight), we will set them to zero in the investigations to 
come. 

As for second order transformations, determined by the matrix A, they can be 
conveniently mapped into the group Sp2n,n of real symplectic transformations, acting linearly 
in phase space (as second order transformations acting on the Hilbert space make up the 
multi-valued metaplectic representation of the symplectic group |[T9l ). Recall that a matrix 
S belongs to the symplectic group Sp2n,R. if an d only if it preserves the antisymmetric form 
0: S € SX(2n,lR) : S G Sp 2n ,iR, S T ttS = Q. Let us also recall that a symplectic 
transformation S acts by congruence on a covariance matrix a: a i— > S T crS. Of course, 
symplectic transformations can in general affect both the entanglement and the energy of a 
state. The algebra of generators of the symplectic group is comprised of all the matrices that 
can be written as ft J, where J is some 2n x 2n symmetric matrix l20l (in this notation, 
generators are not complexified, so that S = c nj ). Such generators do not have a definite 
symmetry (i.e., they are not necessarily symmetric or antisymmetric). Choosing a basis 
of the algebra such that each generator of the basis is either symmetric or antisymmetric, 
allows one to distinguish between a compact subgroup K(n) = Sp2n,n H SO(2n) (spawned 
by antisymmetric generators) and a non-compact subgroup (arising from skew-symmetric 
generators). Notice also that, since compact transformations are indeed orthogonal, they do 
not affect the energy of the states they act upon (explicitly, Tr cr and ||i?|| 2 are both invariant 
under phase space "rotations"). 

Remarkably, the subgroup K(n) is isomorphic to U(n). Because this fact will be 
exploited throughout the whole work, we shall sketch its proof here. Let us define the 
transformation O by 

X Y 
W Z 



o = 
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where X, Y, W and Z are n x n real matrices. It is straightforward to show that this 
transformation is symplectic and orthogonal if and only if Z = X, W = —Y, X T X+Y T Y = 
1 and X T Y - Y T X = 0, so that 

°=(-y Y x) < 4 > 

Now, let U = X + iY be a matrix with real part X and imaginary part Y. The unitarity 
condition on U corresponds exactly to the previous two conditions on X and Y, thus 
demonstrating the existence of a bijective mapping from U(n) to K(n). The preservation 
of the composition rule can be straightforwardly checked out. Incidentally, this isomorphism 
implies that K(n) has n 2 independent parameters. 

Let us rephrase Eq. (f2]i to give a transparent parametrisation of pure Gaussian states 
in phase space terms, by considering the action of first and second order operations on the 
covariance matrix and on the first moments: 

a = S T S , with S € Sp 2n ,M , R € R 2 ™ (5) 

(recall that, in our units, the covariance matrix of the vacuum is the identity). Indeed, 
because of the peculiar nature of their characteristic functions, Gaussian states are completely 
determined by first and second moments of the canonical operators. 

More generally, let us also recall that the CM X of any, pure or mixed, Gaussian state 
can be written as 

X = S J vS , (6) 

where S G Sp2n,n an d v = diag(Vi, . . . , v n , v\, . . . , v n ) is a diagonal matrix with double- 
valued eigenvalues called the "Williamson normal form" of £ ||2T1 l22l (corresponding to 
the normal-modes decomposition of positive definite quadratic Hamiltonians). The real 
quantities {i/j} are referred to as the 'symplectic eigenvalues' of S and can be computed 
as the eigenvalues of the matrix |ifffi|. The symplectic eigenvalues hold all the information 
about the entropic quantities of the Gaussian state in question. In particular, the 'purity' 
fj, = Tr q 2 of the Gaussian state q with CM S is determined as 



vn 



V 3 
1 



1/VDetS,, (7) 

while the von Neumann entropy S = — Tr (g In g) reads 

n 

s = J2 h ^> < 8 > 

with the 'entropic function' h(x) given by 

K x ) = !og2(— ^— ) — log 2 (^— )• (9) 

Being the eigenvalues of |if2S|, the symplectic eigenvalues are continuously determined by 
'symplectic invariants' (i.e. by quantities depending on the entries of the CM invariant under 
symplectic transformations), defined as the coefficients of the characteristic polynomial of 
such a matrix [24 25 1). This observation will be useful later on. 
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3. Measures on the set of pure Gaussian states 

The present section is devoted to the definition of consistent measures on the set 
of pure Gaussian states, introducing a broad framework motivated by fundamental 
statistical arguments. We will review the construction of the 'micro-canonical measure', 
already introduced in Ref. ifTTl . complementing such earlier studies with discussions and 
mathematical details. Furthermore, within the same general framework, we will present a 
novel, "canonical" measure, thus extending our previous treatment. 

A 'natural' measure to pick would be one invariant under the action of the operations 
which generate the set of states we are focusing on. In the previous section, we have analysed 
such a set of operations for pure Gaussian states, showing that it amounts to symplectic 
operations and displacements. One would thus be tempted to adopt the left- and right- 
invariant measure {i.e., the Haar measure) over such groups. Unfortunately, because the 
symplectic group is non compact, the existence of a Haar measure on the whole group [from 
which a measure for pure Gaussian states could be derived via Eq. (|5]l] is not guaranteed. 
Notably, even if such a measure could be constructed, it would not be normalisable, giving 
rise to distributions with unbounded statistical moments. Moreover, some prescription has 
obviously to be introduced also to handle the first moments which are, in general, free to 
vary in the non-compact R 2n (notice that the Euclidean volume is obviously invariant under 
left translations but is not a proper measure in the space of first moments because it is not 
normalisable, due to the non-compactness of R 2n ). 

To cope with such difficulties we will introduce, in analogy with statistical mechanical 
treatments, assumptions on the energy of the states under examination, which will constitute 
our 'privileged' physical observable (in a sense to be elucidated in the following). A proper 
structure to introduce a measure is inspired by a well known decomposition of an arbitrary 
symplectic transformation S: 

S = O'ZO, (10) 

where 0,0' E K(n) = Sp(2n,lR) n SO(2n) are orthogonal symplectic transformations, 
while 

Z = Z'S)Z'~ 1 , (11) 

where Z' is a diagonal matrix with eigenvalues zj > 1 V j. The set of such Z's forms a 
non-compact subgroup of Sp2 n ,m (corresponding to local squeezings), which will be denoted 
by Z(n). The virtue of such a decomposition, known as "Euler" (or "Bloch-Messiah" l23l ) 
decomposition, is immediately apparent, as it allows one to distinguish between the degrees of 
freedom of the compact subgroup (essentially 'angles', ranging from to 2ir, which moreover 
do not affect the energy) and the degrees of freedom z/s with non-compact domain. In 
particular, applying Euler decomposition to Eq. (f5]l leads to 

<t = T Z 2 . (12) 

Due to the rotational invariance of the vacuum in phase space, the number of free parameters 
of a pure Gaussian state of an n-mode system is thus n 2 + 3n (taking the 2?? independent first 
moments into account). 

Quite naturally, we shall assume the n 2 parameters of the transformation O to be 
distributed according to the Haar measure of the compact group K(n), carried over from U(n) 
through the isomorphism described by Eq. (3). The set of such parameters will be compactly 
referred to as while the corresponding Haar measure will be denoted by dpni^)- 

We have thus identified the variables which parametrise an arbitrary pure Gaussian state 
and imposed a distribution on a subset of them. A 'natural' measure has yet to emerge 
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for the non-compact variables {zj} and {Rj}. In order to further constrain the choice of 
measures on such variables we shall invoke now a fundamental statistical argument. In their 
kinematical approach to statistical mechanics ||26l , Popescu, Short and Winter introduced a 
general principle, which they refer to as general canonical principle, stating that 

"Given a sufficiently small subsystem of the universe, almost every pure state of the 
universe is such that the subsystem is approximately in the 'canonical state' q c ." 

The 'canonical state' g c is, in our case, the local reduction of the global state picked 
from a distribution of states with maximal entropy under the constraint of a maximal total 
energy E. That is, quite simply, a "thermal state", which is a Gaussian state with null first 
moments and CM er c = (1 + T/2)l. Here the 'temperature' T is defined by passage to 
the "thermodynamical limit", that is for n — > oo and E — ► oo, (E — 2n)/n — ► T (assuming 
ks = 1 for the Boltzmann constant). For ease of notation, in the following, the symbol ~ will 
imply that the equality holds in the thermodynamical limit, e.g.: (E — In) In ~ T. Notice 
that we have required here the introduction of a maximum energy E or, alternatively, of a 
temperature T. In point of fact, such requirements are necessary to handle the non-compact 
part of the symplectic group. Note also that, in principle, two options are open in this respect, 
as one can introduce either an upper bound to the energy or a temperature: essentially, these 
two distinct options characterise the two distinct approaches (micro-canonical and canonical) 
which we will detail in the next section. 

Because the canonical state g c is Gaussian with vanishing first moments, the general 
canonical principle can be fully incorporated into our restricted (Gaussian) setting. As we 
have shown in Ref. ifTTI . the compliance with the general canonical principle enforces a rather 
stringent restriction on the distribution of the non-compact variables. In particular, the general 
canonical principle is always satisfied if, in the thermodynamical limit, such variables are 
independent and identically distributed (i.i.d) E7\ . To keep our exposition lighter and more 
readable, we will not repeat here the technical derivation of this implication. Actually, it will 
be entirely subsumed, a posteriori, by the derivation of the "concentration of measure" in 
Sec.g] 

Before moving on with the definition of specific measures, let us comment on the first 
moments {R}. The general canonical principle imposes that, in the thermodynamical limit, 
their density of probability p'(R) tend to a ^-distribution centered in (in fact, the canonical 
state has vanishing first moments). A suitable example is p'(R) = (n\/ir) n exp(— nA|ji?|| 2 ) 
for some constant A (notice that ||i?|| 2 is the first moments' contribution to the energy). Let 
us remark that this class of distributions encompasses the ones usually adopted for coherent 
states, in the computation of classical teleportation thresholds [28 ] 29 j . From now on, we will 
just set the first moments to zero: they can be coherently incorporated into our general picture 
following the recipe given above. 

Let us now turn to second moments and sum up our line of thought so far: inspired by 
mathematical considerations and guided by physical arguments, we have defined a distribution 
for the "compact" degrees of freedom *& (essentially, their Haar measure) and specified a 
prescription for the non-compact parameters z/s. Several choices are then possible, within 
this prescription, to deal with the variables z/s. We will describe now in detail two of such 
choices, which we will then apply to study the typical entanglement. 

3.1. Micro-canonical measure 

As a first approach, we will introduce a micro-canonical measure on the class of n-mode 
pure Gaussian states with an energy upper bounded by E. Notice that such a restriction is 
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essentially equivalent to fixing the total energy of (n + l)-mode states to E + 2 (in fact, in the 
latter instance, the energy of the additional mode is not independent and merely 'makes up' 
to reach the fixed total amount). Notably, the two approaches are obviously indistinguishable 
in the thermodynamical limit. 

Notice that the parameters Zj, whose distribution is left to define, determine the energy 
Ej pertaining to each decoupled mode j of the Euler decomposition, according to (see 
Eq. (O) 

E j = z] + ±. (13) 

z j 

Here, we will assume a Lebesgue ('flat') measure for the local energies Ej's, (uniquely 
determining the squeezings Zj's, as Zj > 1), inside the region Te = {E : \E\ < E} 
bounded by the linear hypersurface of total energy E (here, E = (E\, . . . , E n ) denotes the 
vector of energies, with all positive entries, while \E\ = 2j=i Ej)- More explicitly, denoting 
by dp m c(E) the probability of the occurrence of the energies E, one has 

dp mc (E) = Nd n E = NdE 1 ...dE n if EeT E , 

dp mc (E) = otherwise , (14) 

where Af is a normalisation constant equal to the inverse of the volume of Te ■ Notice that such 
a fiat distribution is the one maximising the entropy in the knowledge of the local energies of 
the decoupled modes. In this specific sense such variables have been privileged, on the basis of 
both mathematical (the Euler decomposition) and physical (analogy with the micro-canonical 
ensemble) grounds. Let us also mention that, as will become apparent in the next section, 
employing the variables {Ej} leads to a remarkable simplification of the expression of the 
averages over the Haar measure of the compact subgroup. While purely formal, this aspect 
yields some significant insight into the privileged role of such variables in characterising the 
statistical properties of physical quantities. 

The micro-canonical average Q mc (E) over pure Gaussian states at maximal energy E 
of the quantity Q(E, i9) determined by the second moments alone will thus be defined as 

Qmc(E) = AfJdp H {d)J^ dEQ(E,#) , (15) 

where the integration over the Haar measure is understood to be carried out over the whole 
compact domain of the variables i9. More explicitly, the integral over the energies E can be 
recast as 

,£-2(n-i) /•E-Ep 1 1 B i 
dE = AT / dE l ... dE n (16) 



T E 

(each energy is lower bounded by the vacuum energy, equal to 2 in the convention adopted 
here), determining the normalisation as Af = nl/(E — 2n) n and leading to a marginal density 
of probability P n (Ej , E) for each of the energies Ej given by 

n ( E i -2\ n - 1 
P n (Ei,E) = 1 3 - . (17) 

1 J ' ' E-2n V E-2n) 

Clearly, the energies Ej are not i.i.d. for finite n. However, as is apparent from Eqs. ( fT6b 
and dTD i, in the thermodynamical limit the upper integration extremum diverges for each 

Ej-2 

Ej while, for the marginal probability distribution, one has P n (Ej,E) ~ c t /T. in 
the thermodynamical limit, the decoupled energies are distributed according to independent 
Boltzmann distributions, with the parameter T playing the role of a temperature, in 
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compliance with the equipartition theorem and equivalence of statistical ensembles of 
classical thermodynamics. This argument shows that the micro-canonical measure fulfills the 
general canonical principle. Also, it naturally brings us to introduce a 'canonical' measure on 
the set of pure Gaussian states. 



3.2. Canonical measure 

In the 'canonical' approach, we will assume for the energies E a probability distribution 
dp c (E) reading 

-(\E\-2n)/T « / -(S 3 -2)/T \ 

dp c (E) = — dE = [] dEj , (18) 

introducing a 'temperature' T. This distribution maximises the entropy on the knowledge of 
the continuous variables Ej's for given average total energy E av , such that riT = E av (the 
latter relation is easily derived by applying Lagrange multipliers). 

The 'canonical' average Q C (T) over pure Gaussian states at temperature T of the 
quantity Q(E, #) determined by the second moments alone will thus be defined as 

r f e ~{\E\-2n)/T 

Q C (T) = I dustf) I — dEQ(E,0) , (19) 

where the integration over the energies is understood to be carried out over the whole allowed 
domain (Ej > 2 V j). 

As already elucidated, the micro-canonical and canonical approaches coincide in the 
thermodynamical limit, as one should expect in analogy with the indistinguishability of the 
classical statistical ensembles in the thermodynamical limit. 



4. Concentration of measure 



In the present section, we will study the statistical properties of the entanglement of pure 
Gaussian states in the thermodynamical limit under the measures introduced in the previous 
section. More specifically, we shall focus on the behaviour of the entanglement of a subsystem 
of m modes (as quantified by the von Neumann entropy of the reduction describing such a 
subsystem), keeping m fixed and letting the total number of modes ?i ^ oo. As we have seen, 
the two measures coincide in this limit (when, in the micro-canonical treatment, the energy - 
notably the other extensive quantity in play - diverges as well). It will thus suffice to consider 
the canonical measure, and the micro-canonical averages will be retrieved upon identifying 
T = (E — 2n)/n. In this section, the shorthand notation x will stand for the average of the 
quantity x with respect to the canonical state-space measure. Under the previous assumptions, 
we will determine the average asymptotic entanglement and rigorously prove that the variance 
of the entanglement tends to zero in the thermodynamical limit. The latter property, which 
had been previewed in Ref. ifTTl . will be referred to as "concentration of measure". 

Let us first recall that the von Neumann entropy S of the 7«-mode reduction is determined 
by the local symplectic eigenvalues {vj, fori < j < m} of the reduced m-mode CM 7 
according to: 

m 

S = J2 h ("o), (20) 

i=i 
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where the entropic function h(x) is defined by Eq. (0, with the additional proviso h(l) = 
(which renders the function continuous). Recall also that the uncertainty principle reads 
i/j > 1 in terms of the symplectic eigenvalues. 

In turn, the symplectic eigenvalues are continuously determined by m symplectic 
invariants {A™}, given by the even order coefficients of the characteristic polynomial of 
the matrix fier J24). Let us denote by g the continuous function connecting the invariants to 
the von Neumann entropy: 

S = 5 (Ar,...,A™). (21) 

Note also that any symplectic invariant A™ is a homogeneous polynomial of order 2d in the 
entries of 7. 

In order to prove the concentration of measure for the entanglement, we will show that 
the distribution induced by the canonical measure on the space of the symplectic invariants 
tends, in the thermodynamical limit, to a 5-function centred on their averages (which will also 
be determined). Through Eq. (ETT i. this will allow us to infer concentration of measure for the 
von Neumann entropy (and to determine its asymptotic average as well). Notice that this will 
also imply concentration of measure for any other entropic measure as, for Gaussian states, 
all such quantities are univocally determined by the symplectic invariants l24l [30l . 

Let us consider the CM <r of the whole system. Considering the expression (fT2l for 
a pure covariance matrix and parametrising the transformation O E K(n) in terms of the 
matrices X and Y such that (X + iY) = U € U (n), according to the isomorphism of Eq. (0J, 
one has for the entries of cr: 

a jk = (XjlXklz? + Yjtfkizf 2 ) for 1 < j, k < n , (22) 
1 

<Jjk = ( Y 3i Y ki z ? + XjiXkizf 2 ) for n + 1 < k < 2n , (23) 
I 

(Tjk = {- X 3i Y kizt + YjiXkiz" 2 ) for 1 < j, (k - n) < n . (24) 

As previously remarked, each of the symplectic invariants is a homogeneous polynomial in 
such entries. Let us now consider, for fixed squeezings Zj's, the averages of such polynomials 
with respect to the matrices X and Y, distributed according to the Haar measure over U (n). 

To work out averages over the Haar measure we shall make use of some basic properties 
of the integration over the unitary group, derived from simple symmetry arguments. The 
reader is deferred to Ref. BTl for a general discussion of such strategies. Let us also mention 
that an alternative way for computing Haar integrals with applications to linear optical systems 
and average entanglement has been discussed in Ref. 0321 . In particular, since permutation of 
the indices is clearly a unitary operation, and since the Haar measure is both left- and right- 
invariant with respect to any unitary operation, the integration over the group only depends 
on the number and multiplicity of the different left and right indices present, but not on their 
specific values. Likewise, the measure is invariant under local phase changes (i.e. operations 
represented by diagonal matrices with one eigenvalue equal to e 1 ^ and all the other eigenvalues 
equal to 1), which implies that all the Xjk's can be swapped with Y^'s without affecting the 
values of the integrals. Another consequence of the invariance under local phase changes will 
be exploited shortly. 

Let us consider the average A™ of the invariant A™ in the limit n — > 00. This invariant 
is a polynomial of order 2d of elements given by Eqs. d22b-d24l), so that each term of the 



Canonical and micro -canonical typical entanglement of continuous variable systems 



12 



polynomial is a product of sums, of the form 

2d 2d / n \ 

n^=n E^L < 25 > 

s=l s=l / 

where Wjki are polynomials of order two in X and Y depending on the z;'s and determined 
by Eqs. (f22b-(f24l>. Suppose we expand the product occurring on the LHS of Eq. ( 1251 ). and sort 
the contribution of the resulting addenda according to the multiplicity of the different indexes 
I's, for 1 < I < n. As already remarked, the integral over the Haar measure of the monomials 
in the matrix elements Xji and Yji, depends only on such multiplicities. Elementary 
combinatorial arguments show that the contributions of leading order in n are the ones for 
which all the indexes Ik's are different from each other, which encompasses n\/(n — 2d)! 
terms of the resulting sum: all the other terms can be neglected in the thermodynamical limit. 

Indeed, only one of the leading terms is non-vanishing. Whenever a term of the sum of 
Eq. ( f24l > enters in a monomial where all the Ik's are different, the resulting average is of the 
form XjiYk\q(X, Y), where the index 1 has been fixed with no loss of generality (due to 
permutational invariance), while the q(X, Y) is any function of the entries of X and Y where 
Xji and Yfci do not appear. In terms of the complex unitary matrix U — X + iY, one then 
has 

teXjtYnq&Y) = (U^ + U^iU^-U^qiX^Y) = , 

because the average has to be invariant under arbitrary left and right phase changes on the 
indexes j, k and 1, but all the 'phase-invariant' terms \Uij | 2 get cancelled out [this holds also 
for j = k, as in that case the factor depending on Uji reduces to (Uj 1 — U* 2 )]. Moreover, if 
all the Ik's are different, terms coming from the sums of Eq. (1221 are identical to those coming 
from the sums of Eq. (l23l l. as X's and Y's indexes can be swapped. Summing up, the leading 
term as n — > oo of the average A™, reads 



A™ ~ c(d, m,n)J2 II ^ + ^ = c ( d ' m > n ) X) II El ' (26) 

where the sum runs over all the possible 2d-subsets of the first m natural integers {i.e. over 
all the possible nl/(n — 2d) ! ~ n 2d combinations of 2c? integers smaller or equal than n, with 
no repetitions, equal to all the possible combinations of different indexes), c(d, m, n) is a 
factor depending on d, m and n which we shall determine shortly, and the symbol ~ specifies 
that the equality holds in the thermodynamical limit. Note that, being interested in finite 
subsystems in the thermodynamical limit, we can always assume 2d < n. The convenience 
of the parametrisation of the states through the energies {E^} becomes fully apparent in 
the previous expression. The coefficient c(d, m, n) is easily determined by considering that, 
setting Ek = 2 Vfc, one has er = 1 regardless of the applied orthogonal transformation, for 
which the values of the local invariants A™ are fixed and trivially equal to their averages. This 
yields 



A™ ~ zM 2d (71) 



' m \ 52s% d Tiles?, E i 
d) ^ 

Notice that the factor ■nT 2d in the normalisation is due to the integration over the Haar 
measure. In fact, such a dependence is common to the average of any monomial of order 
2d in the matrix elements (when non null) ||3T1 . Also, let us emphasise that the previous line 
of thought applies to any polynomial in the entries of <x. The only specific information about 
the invariants A™ entered in imposing the normalisation condition, so that all the structure of 
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the invariants could be 'extracted' from a trivial situation (one has A™ = (™) for the vacuum 

123). 



The same arguments apply to the average (A™) 2 (as a squared invariant is just a 
polynomial of order Ad in the entries of er), leading to 



m\ Y,s™ Ihe s™ E i 
d) ' (2n) 4 



To prove the concentration of measure for the invariants, we have to average over the 
remaining variables {Ej}. As we have seen, both the canonical and the micro-canonical 
distributions are i.i.d. in the thermodynamical limit. For the sake of clarity, let us impose the 
two conditions of invariance under permutations and statistical independence one at a time. 
The invariance under permutations of the variables {Ek} allows one to write Eqs. d27l ) and 
as 



A' 



■in 



2d)\(2n) 2d ~ [d 2 2d V 



d " \dj (n - 2d)\(2n) 2d \d J 2 2d 



n 



leading to 



i ( a ™ 2 - ( a ™) 2 ) = % l ( n e » ) - n ^ I < 29 > 




Notice that the quenching of the residual terms in this development is due to the integration 
over the Haar measure of the compact symplectic group, which always results in terms of the 
order n~ 2j for polynomials of order 2j in X and Y. Finally, the statistical independence of 

the energies Ej's requires n?=i Ej = E k V/c, such that the previous equation becomes 



lim (A" 1 - (A" 1 ) 2 ) = . (30) 

n— *oo 

Therefore, the variance of any local symplectic invariant A™ vanishes in the thermodynamical 
limit. The previous argument applies to any polynomial of finite order in the entries of a. In 
particular, in the case 2d = 1, this vanishing of the variance directly applies to the entries 
of er, thus also implying, a posteriori, the compliance (already shown in Ref. ifTTl in a more 
specific background) of the presented measures with the general canonical principle. 

To complete our proof and extend the concentration of measure to the local von Neumann 
entropy, let us consider the measure r„(A™, . . . , A™) induced, in the m-dimensional real 
space of the local symplectic invariants {A™}, by the proposed canonical measure on n- 
mode pure Gaussian states. Because of Eq. (l30l l. such a measure completely concentrates in 
the average values A™. Here, we will express this fact by claiming that r„ tends to a Dirac 
delta: 

lim r„(Ar, . . . , A™) = S(A[ n - Af , . . . , A™ - Af ) . 



[see Appendix B for a rigorous formulation, showing the emergence of the following limits 
from Eq. d30bl - For the von Neumann entropy of the m-mode subsystem S m , one finds 

lim = lim f r„(Af, . . . , A-)g(Ar, . .., A™) cTA^ = ff (Sf , . . . , A™) , (31) 



lim SI - Sj = . (32) 
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The asymptotic average S m can be determined as well. First, let us determine the 
canonical averages A™ of the invariants according to Eqs. d27l ) and dTHb . finding 



1 d 

n — >oo 



ra\ 1 



IimA3*=( + . (33) 



Next, let us recall the expression of the symplectic invariants in terms of the symplectic 
eigenvalues {i/j} J24): 

a ™=e n v i- w 

Straightforward substitution in Eq. (l34l shows that the values Vk = (1 + T/2), VI < k < m, 
account for the asymptotic values of the invariants given by Eq. ( f33b . Thus, for the symplectic 
eigenvalues {vj}, one has 

T 

lim 77^ = 1 + — V 1 < fc < m . (35) 

n — >oo 2 

In the thermodynamic al limit, both the introduced distributions, canonical and micro- 
canonical, comply with a form of equipartition theorem: an average 'energy' equal to T/2 
is allotted to each decoupled quadratic degree of freedom. Finally, one has 

— T 

lim S = mh(l + -) . (36) 

n — >oo 2 



5. Study of the typical entanglement of pure Gaussian states 

In the present section, we present a detailed study of the typical bipartite entanglement for 
Gaussian states with a finite number n of total modes. More specifically, we will study 
the statistical properties of the entropies of a reduced subsystem of m modes, when the 
global states are distributed according to the micro-canonical and canonical measures. As 
we have already mentioned, the quantity Dct cr is a proper quantifier of the mixedness of 
Gaussian states, and thus of their entanglement as well, when computed for the reduced 
subsystem of a pure state. For such an entropic quantity, we were able to derive analytical 
expressions for the canonical and micro-canonical average and standard deviation, which will 
be presented in subsection 15. II In subsection 15. 21 we describe a strategy - based on a simple 
linear minimisation - to derive, from such analytical results, exact information about the 
micro-canonical statistical properties of the actual entropy of entanglement (defined as the 
von Neumann entropy of the subsystem). In subsection !5.3l we complement this analysis with 
numerical results, also covering m-mode subsystems. 



5.7. Purity of single-mode subsystems 

For a single-mode subsystem (corresponding to m = 1 in the previous section's notation), 
we have worked out analytically the average and variance of the inverse squared purity 
fi~ 2 = 1/Tr [q 2 ] 2 , which coincides, for Gaussian states, with the determinant of the reduced 
covariance matrix 7. Let us mention that, in this particular instance, such a measure is a 
perfectly legitimate entanglement quantifier, as it induces, on the set of single-mode states, 
the same hierarchy as the von Neumann entropy (in fact, both quantities are determined by 
the single symplectic eigenvalue, see Section[2]i f34l . 
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Making use of the methods described in Ref. BTI , we determined the first and second 
statistical moments over the Haar measure (respectively denoted by n^ 2 and /J^ 4 ) of the 
determinant Dct 7 = /% 2 , to find 

- E 3 E k 2 



-2 



^ 4(n + l)n n+l' ; 



(37) 



_ 4 _ 1 (n-1)! 



16 (n + 3)! 



53 EjE k EiE m + 8 ^ E* E k E t + 12 J2E]E t 



(96 + 16(n - 2)) 53 £j-Efc - 32(n - 1) 53 ^1 + 12Sn ( n - l ) + 384n 



(38) 



where the averages over the variables {Ej} are still to be worked out, according to the chosen 
distribution. 

The canonical and micro-canonical averages are then straightforward to compute. In the 
canonical instance, the average and second moment u~ 2 and u~ 4 read: 



,- = 1^ + 4^ + 1 



(39) 



_ 4 = 1 nl{n- 1) 
: ~ 16 (n + 3)! 

+32(n + 3)(n + 2)T 



{n 2 + IIti + 22)T 4 + 8(?i 2 + 8n + 6)T 3 + 8(3n 2 + 15n + 10)T 2 

1. (40) 



Whereas, in the micro-canonical case, defining E = (E — 2n), one gets 
(*»-!) ^2 



-2 



4(?i + 2)(n + l) 2 



£ 2 + 4(n + 2)£ +1 



(41) 



-4 



(n!) 2 (n-l) 



((ti 2 + lln + 22)£ 4 + 8(71 + 6)(n + 4)(n + 1)£ 3 



16(n + 4)!(n + 3)! 

+8(n + 4) (71 + 3)(3n 2 + 15n + 10)£ 2 + 32(n + 4) (71 + 3) 2 (n + 2) 2 E ) + 1. (42) 



Notice that the maximal u M for given energy E = E + 2n is (see Appendix C 1 
{E + A) 2 



Mm 



16 



(43) 



Whilst restricted to a particular quantifier and to a single mode, these analytical results 
display the most relevant statistical features of the entanglement of pure Gaussian states. The 
micro-canonical mean /j^ 2 . is monotonically increasing with E for fixed n and, for n > 2, 
monotonically decreasing with n for given E (a "finite size" effect shows up for E < 10, 
where /i^ 2 increases in going from 2 to 3 modes). This can be promptly explained as more 
available energy generally allows for higher entanglement, while the presence of more modes 
"drains" energy away to establish correlations which do not involve the particular chosen 
mode. On the other hand, the canonical average entanglement is monotonically increasing in 
both the temperature and the number of modes. This behaviour is encountered also for the 
micro-canonical entanglement with given maximal total energy per mode, which is, even for 
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small n, closely akin to the canonical ensemble (upon replacing E/n with T). The increase of 
the average canonical entanglement with increasing number of modes but fixed temperature 
is a non trivial, purely 'geometric' effect, due to the average over the Haar measure d/i#($) 
of the compact variables An analogous increase is in fact observed assuming a given, fixed 
value for the variables Zj's and averaging only over the compact variables as the number of 
total modes increases, a given mode has more possibilities of getting entangled, even keeping 
a fixed mean energy per mode. 

As for the standard deviations, which are straightforward to derive from the expressions 
above, they are generally increasing with total energy and temperature for fixed total number 
of modes (as more energy allows for a broader range of entanglement). Significantly, these 
partial analytical results clearly show the arising of the concentration of measure around 
a thermal average. Both for the canonical case and for the micro-canonical one with 
E = nT the standard deviation decreases with increasing number of modes n, falling to 
zero asymptotically (after transient "finite size" effects, for very small n). Moreover, in the 
micro-canonical instance, the thermal average of concentration is generally very distant, even 
for relatively small n, from the allowed maximum of Eq. ( l43b (which clearly diverges in the 
thermodynamical limit): e.g., for E = lOn one has that the average /x~^ is, respectively 16.5 
and 257.1 standard deviations away from the maximal value p^f for n — 5 and n = 20. Such 
a distance increases monotonically with the total number of modes. This peaked concentration 
for finite n will be exploited in the next subsection to obtain strict bounds on the average von 
Neumann entropy of entanglement of single-mode subsystems. 



5.2. Estimating the micro-canonical mean entropy of entanglement 

In the previous section, we have focused on the inverse square purity p~ 2 , because it can be 
readily described in terms of the covariance matrix's entries, thus allowing one to determine 
analytical expressions for the entanglement's statistics. Let us stress once more that, for 
pure states, such a quantity is a legitimate entanglement measure, as it is a monotone strictly 
related to the so-called linear entropy, given by 1 — p. Moreover, for single-mode Gaussian 
states, the linear entropy is, as we have already remarked, monotonically determined by the 
von Neumann entropy alone. However, the proper "entropy of entanglement" (given by the 
von Neumann entropy of the reduced density matrix) is endowed with a clear operational 
meaning for pure states (as it corresponds to the rate of singlets distillable from the state in 
the asymptotic limit of infinite copies [33]). It is thus highly desirable to obtain qualitative and 
quantitative information about the statistical properties of such a quantity under our measures. 
Clearly, the findings of Section[4]already provide us with such results in the thermodynamical 
limit, where the von Neumann entropy concentrates, with vanishing variance, around the value 
set out in Eq. d36b . 

Here, we shall address the entropy of entanglement of a single-mode subsystem in a 
system with a finite total number of modes. Exploiting the fact that the von Neumann entropy 
S is continuously determined by p~ 2 as l34l 

S(p) = hip- 1 ) (44) 

(under the additional prescription 5(1) = 0), we will show how upper and lower bounds on 
the micro-canonical average S mc can be derived. The reason why we focus on the micro- 
canonical measure is that, imposing a bound on the energy, it involves a maximal value for 



the entanglement (as discussed in Appendix C I, which is another key ingredient we are going 



to use. Extending the method we are presenting to the canonical measure is possible, at the 
price of introducing an approximation to neglect the 'tail' of the canonical distribution in 



Canonical and micro -canonical typical entanglement of continuous variable systems 



17 



1 


3 


1 


2 


1 


1 




1 





9 





8 





7 





6 



♦ j * 

♦ ★ 


♦ ★ 






♦ * 


Entropy of entanglement 




★ 


♦ upper bound 




♦ 


★ lower bound 




* 







6 8 10 12 14 16 18 20 
n 



Figure 1. Upper and lower bounds to the micro-canonical average von Neumann entropy of 
entanglement of one mode as a function of the number of total modes n. The maximal energy 
was assumed to be An (in units of tkj/4). To compute the bounds, the method described in 
subsection !5.2l has been employed. The parameter M was chosen to be 10000. 



the space of the entanglement (certainly feasible with reasonable resources for small enough 
temperatures). As we will see, such bounds, which we know to become increasingly close 
with increasing n because of concentration, can be quite tight even for relatively small n, thus 
providing precise information about the average entropy of entanglement. 

For the sake of readability, we will from now on set a := /i~ 2 . For a given number 
n of modes and total energy E in the system, the maximum value a max for the inverse 
squared purity of any fixed mode is given by Eq. d43l . while the minimum value is always 
(for n > 2) given by a m in = 1, corresponding to the case where the mode is decoupled 
from the remainder of the system. Our measure on Gaussian states will induce a probability 
distribution v on the interval [a m i n , a max ]- While we do not know v directly, we are aware of 
two of its properties, namely the averages a c and a 2 with respect to it, which we computed in 
the previous subsection. Our approach is going to be as follows: to obtain an upper bound for 
Smc we maximise S mc over all probability distributions v ' which produce the right averages 
for a and a 2 . There is a technical obstacle to the pursuit of this programme: the set of all 
probability distributions on the interval is infinite-dimensional. To circumvent this problem, 
we will partition [a m ; n , a max ] into M equally sized sub-intervals (M £ N being an arbitrary 
parameter) and consider "discretised" probability distributions, which are constant over these 
subintervals. This is going to be done in such a manner that the resulting bounds are valid for 
any finite M, and not only in the limit of M — ► oo. Also, it will turn out that all optimisations 
can be cast into the form of linear programs. This excludes the occurrence of local minima 
and implies that the obtained bounds will hold rigorously. 

More explicitly, fix an M £ N. The kth sub-interval will be [l(k),l(k + 1)], where 
l(k) = a m j n + k a ™»x^a m in is its leftmost point. Let := v([l(k), l(k + 1)] be the measure 
of the fcth sub-interval. Clearly, one has 

/-a max M-l M-l 

a mc = j aAv{a) < l(k + l)vk, a mc > ^ K k ) v k , (45) 

^ a """ k=0 k=Q 

M-l M-l 

fc=0 fe=0 
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Figure 2. Histograms of microcanonical entanglement distributions, m = 1, and 5000 
samples were taken for each distribution. The concentration of measure with increasing n 
is apparent. 



Furthermore, because S(a) is monotonously increasing in a, we have 



Af-l 



S mc < sup V S(l(k + l)y k . (48) 

,,'IZP '~ d 



k=0 



S mc = f S(a) di/(o) < ^ + 1)H • (47) 

Now let P C R A/ be the set of probability distributions which are constant on the M sub- 
intervals and compatible with Eqs. d45b . d46b . Certainly, the discrete distribution given by the 
Vf. is an element of P and we can thus deduce from Eq. (l47l i that 

Af-l 

v'ep 

Hence, the following linear program with variables v' € R M yields an upper bound for S„ 

maximise S(l(k + l))vk 

k 

subject to a mc < y~]l(k + 1)^, a mc >^Z(fc)^ , 
fe fc 

A.- fc 

fc 

A lower bound is found by a completely analogous minimisation program. 

Figure Q] depicts the results of the programme laid out above for a specific choice of 
energy per mode. Quite remarkably, already for a small number of modes, one can give a fairly 
precise value for the average micro-canonical entropy of entanglement, thus complementing 
the asymptotic statement of Section|4] Let us also mention that the restriction imposed by the 
knowledge of the second moments is crucial in rendering the bounds so tight (as opposed to 
the use of the averages alone). 
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Figure 3. Numerical evidence that the ratio variance/average, for a microcanonical measure 
on n-mode states with maximal energy E, is proportional to 1/n and, surprisingly, that 
the proportionality factor is independent of E/n. A single-mode subsystem was considered 
(m = 1); 1500 states were sampled for each data point. 



5.3. Numerical results on the entropy of entanglement 

Unfortunately, the analytical and exact methods illustrated above cannot be easily extended 
to more complicated situations. However, our measures are well suited to be numerically 
investigated by direct sampling (the Haar measure-distributed symplectic orthogonal can be 
reproduced by generating unitary matrices from the Gaussian Unitary Ensemble l35ll and 
then by translating them according to Eq. This allows for the investigation of the 

statistical properties of the actual entropy of entanglement for varying values of m, n and 
of the parameters of the measures E or T, according to the setting in question. 

In Figure|2]a sequence of numerically generated microcanonical probability distributions 
for the entanglement of a single mode are plotted, unambiguously showing the concentration 
of measure for the von Neumann entropy at small n. Notice that, in the microcanonical case, 
even for small n - well before the onset of thermodynamical concentration of measure around 
the finite thermal average - the entanglement of pure Gaussian states distributes around values 
generally distant from the finite allowed maximum: e.g., for m = 1 and E = 10??., the 
difference between the maximum and the average S is, respectively, 4.0 and 13.6 standard 
deviations for n = 5 and n = 20. 

Furthermore, for the micro-canonical case, our investigation through sampling indicates 
that the ratio between variance and average is inversely proportional to n for finite n and, 
rather surprisingly, that the proportionality factor is independent of E/n (see Fig. |3). 

When more than one mode is addressed (i.e. when m > 1), numerics show that 
the average canonical entanglement is roughly proportional to the number of local modes 
m. The same linear approximation carries over to the micro-canonical case for E 3> n 
and m n, upon substituting E/n for T. Clearly, such an approximate proportionality, 
rigorously true in the thermodynamical limit [see Eq. 1361 1. is better satisfied for increasing 
number of total modes n, but provides good estimates (up to some percent) already for 
n « 30 and small enough m, as is shown in Fig. |4(a)[ where a clear subadditive behaviour 
for m < n is also apparent. Also, let us point out that the concentration of measure for the 
entanglement distribution clearly shows up at finite n for m > 1 as well. For instance, for a 
microcanonical distribution with m = 5, n = 20 and E = 200, one finds that the average 
is 11.4 standard deviations away from the allowed maximum. Also, Fig. |4(b)| shows how the 



| The MATLAB code we made use of is available at www.imperial.ac.uk/quantuminformation. 
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Figure 4. Microcanonical average entropy of entanglement (a) and standard deviation of the 
entropy of entanglement's distribution (b) for E = lOra, as a function of the number of modes 
m of the reduced subsystem. On the left (a), the upper curve (red) refers ton = 30, the middle 
curve (blue) to n = 20 and the lower curve (green) to n = 10. On the right (b), the upper 
curve (green) refers to n = 10, the middle curve (blue) to n = 20 and the lower curve (red) 
ton = 30. Only the range 1 < m < 5 was considered for n = 10. Data were obtained from 
samples of 5000 states. 



standard deviation increases approximately proportionally to m and decreases for increasing 
n for m^las well. 

5.4. Concentration of measure for infinite subsystems 

So far, we have considered the concentration of measure occurring for a finite subsystem 
with a fixed number of modes m. Let us now briefly turn to the case of a fixed 'ratio ' of 
subsystems, where m/n ~ a > in the thermodynamical limit. Clearly, in such a case, the 
subsystem as well is comprised of infinitely many modes 

The analytical reasoning of Section|4]can be readily adapted to this case and shows, not 
surprisingly, that the average S of the von Neumann entropy diverges in the thermodynamical 
limit. In this case though, analogous arguments imply that the variance S 2 — S of the 
von Neumann entropy diverges as well. Still, numerics strongly suggest the occurrence of 
a weaker form of concentration, namely 

(S 2 -S 2 )/S 2 ~0. (49) 

This behaviour, unambiguously supported by the numerical analysis in both the micro- 
canonical and canonical instances, extends to the general case of large m the evidence 
depicted in Fig. [3] (according to which the ratio between variance and average is inversely 
proportional to n). 

6. Conclusions and outlook 

The approach we have introduced here tames the divergence due to the infinite dimension 
of the Hilbert spaces of continuous variable systems and allows one to introduce a well 
defined notion of "typical" continuous variable entanglement. Our exhaustive analytical 
and numerical study shows that such a typical entanglement concentrates sharply around a 
thermal average, even for a relatively small number of modes. In the micro-canonical case, 
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where the upper bound to the energy implies the existence of an absolute upper bound for the 
entanglement, the average entanglement turns out to be different and well separated from the 
maximal value, by many standard deviations, even for very small number of modes. Under 
such 'heavily thermodynamical' prescriptions - as are the ones we have adopted to construct 
the canonical and micro-canonical measures - equipartition prevents the entanglement from 
reaching the allowed maximum, evenly spreading the correlations between all the modes. 

Our measures, being compliant with the general canonical principle, should be suitable to 
the description of dynamical situations l36l|T8l , also involving randomised interactions J7]. A 
systematic study of such processes - notably restricting to two-mode interactions, in the spirit 
of Ref. [7| - where the two measures defined here would arise as stationary distributions, is 
the next direction to pursue in the present line of research. 

This line of enquiry may also contribute to our understanding of why objects typically 
appear classical despite being governed, as is normally assumed, by quantum mechanics. 
In particular, this and other related works, like Refs. J6j [26] [36), imply that, always in a 
particular sense, most global pure quantum states have the property that a local state of a 
comparatively small party, obtained by tracing over the larger party, will be highly mixed and 
accordingly have little or no quantum correlations. This is at least one necessary ingredient 
for the system to appear classical. The work here indicates that these arguments can be 
applied in the continuous variable setting too. It would be interesting to combine, compare and 
further develop existing approaches in order to explain the emergence of apparent classicality 
in realistic models of physical systems. As a further investigation, one could consider the 
extension of analogous micro-canonical and canonical measures on general finite-dimensional 
states [ 18 1. One could then distinguish, to a greater detail, the features induced by the chosen 
thermodynamical setting from the ones proper to the Gaussian continuous variable scenario. 

Finally, let us mention that the presented framework may be suitably extended to more 
general, mixed Gaussian states. In fact, the Williamson decomposition of a generic covariance 
matrix [Eq. ©], together with the Euler decomposition given by Eq. dTOb . suggest that, to 
encompass mixed states as well, one has to add only another set of compact variables d, plus 
the symplectic eigenvalues of the global mixed states {vj}. The presented approach thus also 
paves the way for the definition of more general measures on the whole set of Gaussian states. 



Appendix A. A Baker-Campbell-Hausdorff formula for quadratic polynomials in the 
canonical operators 

For the sake of completeness and self-consistency, we present here a proof of Eq. (0. To 
this end, we will follow a strategy customarily adopted in the derivation of Baker-Campbell- 
Hausdorff-like relations (see, e.g., Ref. 11371 ). Let us thus define the operator GU,&($) = 
e "?(fl AR + R b \ where $ is a real variable. Note that the CCR straightforwardly imply 

[R T AR, R T b] = 4,iAQb . (A.l) 

We start by assuming the following ansatz, which we will prove shortly, 

G A b(tf) = e M{RTAR+RTb ' ) = e ifW(R T A&) e R T (M(#)b) ^ (A 2) 

where /(<!?) is a scalar real function, whereas M($) is a 2n x 2n matrix depending 
continuously on Differentiating both sides with respect to d, we find (differentiation is 
denoted by 
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= (f\0)& r Ail + & T e i 'W&^M\0)be- i 'W&^>) 6 A , b 

= (f'(tf)R T AR + R T c- 4fWAn M'($)b} G A ,b , (A.4) 

where we have made use of Eq. (IA.U in the last step. Equating dA.3b and dA.4t yields the 
following systems of differential equations: 

f'W =1, (A.5) 
M'(tf) = e 4/ ^ An , (A.6) 

with initial conditions /(0) = and M(0) = 0. This system always admits an analytical 
solution, given in general by /($) = i? and M(i?) = J Q e 4l? An di9', whose value in i9 = 1 
gives the matrix M = M(l), thus proving the validity of Eq. ||3). If A is invertible, M is 
simply given by 

M = -^lA- 1 (1 2 „ - e 4AO ) . (A.7) 
Appendix B. Asymptotic concentration of measure 

In this appendix, we show how the concentration of the measure T n in the space of the 
symplectic invariants follows from Eq. (130b . In doing so, we will provide the reader with 
a formal derivation of Eqs.(l3Tb and ( l32l . 

Let us consider the m-dimensional real space A m of the vectors of symplectic invariants 
A = (A™, . . . , A™), endowed with the usual Euclidean norm j| • ||. Let D e be a spherical 
ball of radius e centered in A = (Af , . ._, A™): D e = {A : || A - A|| < e}. Let R e be the 
complement of D e : R e = {A : ||A — A|| > e}. Recall that r n stands for the (normalised) 
measure induced on the space A m by the canonical measure of n-mode pure Gaussian states. 
Also, r n (dA) will stand for the infinitesimal element of such a measure. Let us remark that 
Eq. d30l ), holding V d, implies 

lim (]|AP - || A|| 2 ) = || A - A|| 2 = . (B.l) 

n — >oo 

Our first aim is deriving a rigorous formulation of "concentration of measure", i.e.: 

Ve > and V£ > , 3 n \ Vn > h : T n {R e ) < £ . (B.2) 

To this aim suppose, ad absurdum, that the latter statement did not hold. Then, Vn, 3 e, £ > 
and no > n such that one has T ng (R e ) > £. But this would imply that, Vn, 3 no > n such 

that || A — A|| 2 > £e 2 , which would overtly contradict Eq. ( IB. 11 1. Because of normalisation, 
the following equation, complementary to the previous one, holds as well 

Ve > and V£ > 0, 3 n \ Vn > n : T n (D e ) > (1 - ^) . (B.3) 

Together with Eq. ( IB. lb . the two previous statements entail 

Ve>0andV£>0, 3n|Vn>n: [ \\ A|| 2 r„(dA) < £ . (B.4) 

The consequences of these facts for the local von Neumann entropy can be easily derived 
by exploiting the following simple properties. As apparent from Eq. ( |34l . A™ > {vj 1 ) 2 , 
Vj and V d: any symplectic invariant is larger than any symplectic eigenvalue. Moreover, 
for the function h(x) - defining the von Neumann entropy according to Eq. (l20b - one has 
x 2 > h(x) for x > 1 (as is the case for the symplectic eigenvalues, lower bounded by 1 
because of the uncertainty principle). This results into || A|| 2 > S. One can also show that 
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j| A|| 2 > S 2 . Therefore, bounds analogous to ( |B.4t hold for the integrals over R e of S and 
S 2 as well. Let us also note that the function g(A), relating the symplectic invariants to 
the von Neumann entropy, is certainly continuous (as it relates a continuous function of the 
eigenvalues of a strictly positive matrix to the coefficients of the characteristic polynomial of 
the matrix). Putting everything together we find 

V£ > , 3 n s.t. V?7. > n : (g(A) - - £) + £ < S < (g(A) + £) + £ , (B.5) 
which is just equivalent to Eq. ( f3Tb . In the previous inequalities the continuity of g(A) has 
been invoked and the integral giving the average S has been decomposed into an integral over 
D e and an integral over R £ . An identical argument holds for the average S 2 , which can be 
shown to converge to g(A) 2 , thus proving Eq. (l32l as well and completing our treatment. 

Appendix C. Maximal entanglement for given energy 

We derive here an expression for the maximal value of the entanglement d43l of pure Gaussian 
states for given energy (under a generic m + n mode bipartition), by adopting an explicit phase 
space approach. To begin with, let us remark that any pure Gaussian state of m+n modes can 
be reduced, by local (with respect to the m + n mode bipartition) symplectic operations, into 
the tensor product of m two-mode squeezed states and of n — m uncorrected vacua (here we 
assume, without loss of generality, m < n) |38l[39l . Now, the local reduction of such a state 
pertaining to the m-mode system is a Gaussian state with CM in Williamson form. We will 
now prove that, amongst the CM's with the same symplectic spectrum, the Williamson form 
is the one for which the second moments' contribution to the energy E = Tr er is minimal. To 
this aim, let us recall that a generic CM er with Williamson form v (and symplectic spectrum 
given by {vj, for 1 < j < to} can be written as 

er = O lT ZO T uOZO' , (C.l) 
where the Euler decomposition of a generic symplectic transformation has been applied. 
Clearly, the transformation O' do not affect the energy and can be neglected in what follows. 
As for the other terms, let us define Z' = diag(zi, . . . , z m ), v' = diag(;/i, . . . , v m ) and X 
and Y such that (X + iY) E U (to) [by virtue of the isomorphism of Eq. @], to obtain: 

Trcr^Tr [{Z 12 + Z'- 2 ){X T u'X + Y T v'Y)} > 

> 2Tr [X J v'X + Y T v'Y] = 2TlV = Ttv , 

where the inequality ensues from a basic property of the trace of a product of positive matrices 
(see l40l and notice that, obviously, the eigenvalues of (Z' 2 + Z'~ 2 ) are larger than 2) and 
from the fact that the transformation parametrised by X and Y is orthogonal and preserves the 
trace. The state achieving maximal entanglement for given energy E is thus a tensor product 
of to two-mode squeezed states (being the state with minimal energy for given entanglement). 
The von Neumann entropy S of the rn-mode reduction of such a state is given by Eq. d20l i. with 
the local symplectic eigenvalues {v^ subject to the constraint E = 45^- =1 Vj + 2(n — to). 
Because of the concavity of h{x), the optimal choice of i/j's, maximising the local entropy, is 
simply given by Vj = " E ~ 2 | T "~ ? "" > Vj, in compliace with the previous constraint. Finally, the 
maximal von Neumann entropy S max (m, n, E) of an m-mode reduction of a (to + n)-mode 
pure Gaussian state with m < n and total energy E is 

S max (m,n,E) = mhf— — ^ — —^j . (C.2) 

Note that this expression diverges in the thermodynamical limit. The corresponding minimal 
purity is given by Eq. d43l . 
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